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By reformulating the first law of thermodynamics in the fashion of quantum-mechanical operators 
on the parameter manifold, we propose a universal class of quantum heat engines (QHE) using the 
multi-level quantum system as the working substance. We obtain a general expression of work 
for the thermodynamic cycle with two thermodynamic adiabatic processes, which are implied in 
quantum adiabatic processes. We also classify the conditions for a 3-level QHE to extract positive 
work, which is proved to be looser than that for a 2-level system under certain conditions. As a 
realistic illustration, a 3-level atom system with dark state configuration manipulated by a classical 
radiation field is used to demonstrate our central idea. 
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I. INTRODUCTION 

A usual heat engine operates between two heat baths, 
the temperature difference between which completely de- 
termines the maximum efficiency of the heat engine. Cor- 
respondingly, it is inferred that no power can be extracted 
if the two baths have the same temperature. But things 
become different when we use quantum matter as the 
working substance. Quantum effects can highlight the 
thermodynamic differences between classical and quan- 
tum working substance of a heat engine. Recently, great 
efforts [jj has been devoted to the investigation of the 
quantum effect of the working substance. Some exotic 
phenomena were discovered, most of which concerns the 
following three aspects: The first aspect is whether we 
can improve the efficiency of the QHE to a level be- 
yond the classical limit. For example, M. O. Scully et 
al. 0, 0, Q proposed a quantum-electrodynamic heat 
engine that can exceed the maximum limit by using the 
photon gas as working substance. The second aspect is 
how we can better the work extraction during a Carnot 
cycle. T. D. Kieu present a new type of QHE [5j,|6|, in 
which the contact time is precisely controlled (without 
reaching thermal equilibrium) and the two heat baths 
are specifically modified. This QHE can extract more 
work than the other model in thermal equilibrium case. 
The third aspect is about the constraints of the temper- 
atures of the two heat baths, under which positive work 
can be extracted. It is clear that a classical heat engine 
can extract positive work when and only when Th > 7). 
Here, Th and Ti are the temperatures of the source and 
the sink. But for a QHE, some exoticphenomena may 
occur. One example is given in Ref. 0, in which posi- 
tive work can be extracted even from a single heat bath, 
i.e., Th = Ti\ Another example is the simplest 2-level 
QHE model in Refs. |fj, |?J , in which positive work can 
be extracted only when Th is greater than 2} to a certain 
extent. 
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We consider the third aspect in detail. In Refs. 
the QHE works between the source and sink at tem- 
peratures Th and I) respectively. Within a cycle, the 
level spacing A changes between A h and A'. For such 
a QHE, the system couples to the bath for a sufficiently 
long time until they reach the thermal equilibrium state. 
Then positive work can be performed when and only 
when T h > Ti(A h /A l ) (condition I). This result implies 
a broad validity of the second law and shows by how 
much Th should be greater than 7} such that the positive 
work can be extracted. This constraints about tempera- 
tures is obviously counter-intuitively different from that 
of a classical heat engine. Now we wonder whether it is 
a universal condition for all multi- level QHE. Actually 
T. D. Kieu has considered two special cases, the simple 
harmonic oscillator and the infinite square well, in (the 
appendices of) Ref. 0. As to the two special cases, be- 
cause all the level spacings change in the same ratio, the 
PWC has the same form as that for a 2-level case 0. 
In this article, we will prove that the PWC for a 3-level 
QHE can be looser compared to a 2-level case under our 
criterion when the level spacings change properly in the 
thermodynamic cycle. 



This paper is organized as follows: In section II, we 
formulate a quantum version of the first law of thermody- 
namics for a multi-level quantum system. In section III, 
we discuss the relationship between the quantum adia- 
batic process and the thermodynamic adiabatic process. 
In section IV we analyze the quantum thermodynamic 
cycle of a universal QHE. The obtained results are more 
universal than those obtained from the 2-level system 
II IH EH In section V, we classify the 3-level QHE 
according to the changes of level spacings. We find un- 
der certain condition the PWC for a 3-level QHE can be 
looser than its counterpart for a 2-level case. A realistic 
model — a 3-level atomic system with dark state struc- 
ture manipulated by a classical radiation field — is given 
to demonstrate our central idea in section VI. 



2 



II. QUANTUM VERSION OF THE FIRST LAW 
OF THERMODYNAMICS 

We consider the QHE with a TV-level system as its 
working substance. The system eigenenergies vary adia- 
batically with the parameters R = (i?i j i?2,..., Rm,) in an 
A/-dimensional manifold M. One can manipulate the 
parameters to implement the quantum adiabatic pro- 
cess, which does not excite the transitions among the 
instantaneous eigenstates \n) = |n[R]) of the Hamilto- 
nian H = H\R] with the instantaneous eigenvalues E n 
= \E n [R\). Usually the density operator p is a function 
of both the external parameters R and the thermody- 
namic parameter, i.e., the temperature T. So we need 
to extend the manifold M to include T. For example, 
the heat transfer can be also caused by the change of 
the temperature. In this sense, we define the differential 
one-form on the (M + l)-dimensional "manifold" MT: 
{X = (R,T)} by dF = \7 R FdR+d T F for any function 
F. We emphasize that g?t may be a discrete variation St 
in the time domain. Then in this sense MT is no longer 
a generic manifold. 

During a thermodynamic cycle, work is done by or 
on the system when the mechanical parameters R vary 
slowly. While heat is transferred when the quantum state 
or density operator changes as the temperature varies 
discretely. At two different instants t — and T2, the 
system contacts with the source and sink at temperature 
Th and TJ respectively to reach thermodynamic equilib- 
rium. The couplings of the system to such different heat 
baths yield the discrete change of the probabilities distri- 
bution in every eigenstate. We consider the infinitesimal 
variation 

dU = Tr(pd R H) + Tr(Hdp) (1) 

of the expectation value U = Tr(pH) of the Hamilto- 
nian, where we have introduced the density operator of 
the state p. dU contains two parts corresponding to 
the changes of the Hamiltonian and the density opera- 
tor respectively. Here, the operator 1-form, on the sub- 
manifold M, can be written as 

d R H = d R Q + d R W, (2) 

which is the infinitesimal variation of H. Here, the off- 
diagonal part 

d R Q = {m\dB.H\n)\m)(n\ (3) 
= J2 E »M\n){n\) 

n 

is the heat operator and the diagonal part 

d R W = ^(m\d R H\m)\m)(m\ (4) 

m 

= y^^d R E m )\m)(m\ 

m 



is the work operator. To derive the above equations we 
have used the Feynman-Hellman theorem (m\dH\m) = 
dE m and the formula (m\dH\n) = (E n — E m )(m\d\n) . 

Upon a first glance, the above definitions of the work 
operator d-RW and the heat operator g?rQ are rea- 
sonable. Intuitively speaking, the work can be done 
only when the independent eigenenergy of inner states 
change. This process is just described by c?rW . The off- 
diagonal elements of the infinitesimal variation of H in- 
dicates the transitions among the inner energy levels and 
thus results in the heat transfer. The term d-R,(\n) (n\) 
interprets the heat transferred to or from a quantum sys- 
tem along with the change of the projections to the in- 
stantaneous eigenstate. Physically it means the chang- 
ing of the occupation probabilities rather than the R- 
dependent eigenvalues themselves. In addition, when the 
quantum state p itself changes due to the varying of both 
the control parameters R and the thermodynamic pa- 
rameter T, the heat is also transfered. This kind of heat 
tansfer is described by the second term on the left hand 
side of the Eq. (JTJ. 

In comparison with the existing studies about QHE 6] , 
we write down the infinitesimal energy dU — dQ + dW, 
where the infinitesimal heat transferred dQ and work 
done dW are identified respectively as two path depen- 
dent differentials 

dQ = Tr{Hdp) + Tr(pd R Q) = E ™ d Pm , (5) 

m 

aW = Tr(pd B W)=J2PmdE m . (6) 

m 

Here p m = (m\p\m) are the corresponding occupation 
probabilities in the instantaneous states \m). 

We would like to emphasize that the above two equa- 
tions were even obtained in many previous references 
[H El Ei- But here we present a general derivation, 
the process of which shows the differences between work 
and heat in the view of quantum mechanics. Substan- 
tially, we give the microscopic definitions of work and 
heat based on the non-adiabatic transitions among the 
instantaneous eigenstates of the time-dependent Hamil- 
tonian. 



III. FROM QUANTUM ADIABATIC 
EVOLUTION TO THERMODYNAMIC 
ADIABATIC PROCESS 

The term "adiabatic process" is usually used in both 
quantum mechanics and thermodynamics, and seemingly 
has different meanings for these two cases. Actually, 
a quantum adiabatic process implies a thermodynamic 
adiabatic process. But not all thermodynamic adia- 
batic processes is caused by quantum adiabatic processes 
0, 0| . This understanding is crucial for the following 
analysis about the thermodynamic cycle of the QHE. 



In quantum mechanics, the adiabatic process is de- 
scribed by the time evolution of a quantum system with 
slowly changing parameters. When these parameters 
vary slow enough, the transitions among the instanta- 
neous eigenstates |n[R]) = \n(t)) are forbidden, i.e., the 
system will keep in the n-th instantaneous eigenstate 
\n(t)) if the system is initially in the eigenstate \n(0)) 
at time t = 0. Generally speaking, starting with the 
initial state |^(0)) = ^c„|n(0)), the system will evolve 

n 

into 

W*)> =Vc n expH / E n (t')dt + i ln ]\n{t)}, (7) 
where 

7„ = J\n(t')\^\n(t'))dt' 

is the so-called Berry's phase. This conclusion implies 
that the occupation probability \(n(t) lipit))] 2 in an in- 
stantaneous eigenstate \n(t)} is adiabatically invariant. 
This result is also valid for the initial mixed state, e.g., 
the thermal equilibrium state 

ftQ) = Z>„»K0))(n(0)|, (8) 

n 

where p nn are Gibbs probability distributions. 

We illustrate a quantum adiabatic process and its cor- 
responding thermodynamic adiabatic process in Fig. 1. 
An amount of ideal gas is constrained in a cylinder. The 
piston moves so slow that the gas is always kept in ther- 
mal equilibrium state. Usually no heat is gained or lost in 
a thermodynamic adiabatic process. Thus any isotropic 
process implies its adiabatic property. But all adiabatic 
process are not isentropic. For example, an adiabatic-free 
expansion is not isentropic. Now we prove a quantum 
adiabatic process illustrated in Fig. lb microscopically 
leads to an isentropic process and thus a thermodynamic 
adiabatic process. 
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FIG. 1: (color on line) Illustration of classical adiabatic pro- 
cess (a) and its quantum-mechanical counterpart (6). At time 
t = 0, the piston (wall) locates at 1/(0), and after a long time 
t — t, the piston (wall) moves slowly to L(t). No change 
happens in the energy level populations during the adiabatic 
process. 



invariant during a quantum adiabatic process. Then, 
from Eq. JTUJ, the entropy keeps unchanged. Therefore, 
it is concluded that the quantum adiabatic process of the 
microscopic particles just results in the thermodynamic 
adiabatic process of the macroscopic system. 

The above arguments shows that the quantum adia- 
batic process in quantum mechanics can result in a ther- 
modynamic adiabatic process. But we have to point out 
that not all thermodynamic adiabatic process are caused 
by the quantum mechanical adiabatic process. For exam- 
ple, in the thermodynamic adiabatic process of Fig. la a 
single molecule may experience nonadiabatic transitions 
due to its interaction with other molecules. 



In quantum mechanics, we describe the motion of these 
structureless gas atoms or molecules with an infinite po- 
tential well with one moving boundary (Fig. lb). The 
piston or boundary moves so slow that the quantum adi- 
abatic condition are satisfied ^ij. Now we consider the 
microscopic definition of the thermodynamic entropy 

S= -A;Tr(p[R]lnp[R]). (9) 

The variation of the entropy due to the changes of the 
mechanical parameters R can be calculated as 

lnp„ n + l)rfRP ran . (10) 

n 

According to the quantum adiabatic theorem, the parti- 
cle distributions in the instantaneous energy levesls are 



IV. POSITIVE WORK DONE BY MULTILEVEL 
QHE IN A THERMODYNAMIC CYCLE 

Like the classical heat engine, a universal QHE 
also bases on a thermodynamic cycle. Mathemati- 
cally, it can be understood as a close path I: {X 
H. 7 X 0; = X(t 4 )} on the (A/+l)-dimensional "man- 
ifold" MT. The efficiency of a universal QHE is r\ q = 
AW/Q corresponding to the ratio of the work done to 
the heat absorbed during a cycle, 

aw = jaw (11) 
Q = I 1 (12) 

Jo 
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where t± is the ending point of Step 1. A typical four- 
stroke QHE, consists of two quantum adiabatic and two 
isothermal processes, is a quantum analogue of the clas- 
sical Otto engine, as illustrated in Fig. 2: 




Bath 1, T h Bath 2, T, 



FIG. 2: (color on line) Fhe graphic sketch of a universal QHE 
model. From time t = to t=Ti, the system absorbs heat 
from the source (bath 1). From time t = T\ to t=T2, work is 
done by the system when the piston is pushed by the working 
substance. Step 3 (from time t — T2 to t=T3) and step 4(from 
time t = T3 to t=T4) are reversed processes of step 1 and step 
2 with modification. 



Steps 1 and Step 3 are two isothermal processes. Dur- 
ing Step 1, i.e., from time t = to t = t%, the system 
couples to the source (bath 1) at temperature Th and the 
level spacings keep unchanged. After absorbing energy 
from the heat source, the system reaches thermodynamic 
equilibrium state, which can be described by the density 
operator p(ri). The heat absorbed during Step 1 is 

Q= f 1 J2E m (m\6 T p\m). (13) 

Step 3 is almost an inverse process of Step 1. From time 
i~2 to T3, the system is brought to couple to the sink at 
temperature 7). The density matrix changes into /5(ts) 
after releasing energy into the bath. In the above oper- 
ations, the absorption and the release of energy during 
Steps 1 and Step 3 happen in quantum ways. Namely, 
they can only happen probabilistically. The correspond- 
ing probability depends on the details of the interactions 
and some intrinsic properties, e.g., the temperatures of 
the heat baths. 

Steps 2 and Step 4 are two adiabatic processes. Dur- 
ing Steps 2 the QHE performs positive work when the 
energy spacings decrease. Meanwhile the parameters R 
adiabatically change from R(ti) = R(0) to R(r 2 ). But 
the atomic probability distribution remains unchanged in 



this process. Step 4 is almost an inverse process of step 
2, during which the system is removed from the sink and 
its energy gaps increase as an amount of work is done on 
the system. The net work done by the system during a 
cycle is 

AW = (f ' + f VpA. (14) 

The above results can give the known results in [rl lTolIT^ 
for a 2-level system. 

We assume the thermal equilibrium Gibbs distribu- 
tions for the heat bath. The system will eventually reach 
Gibbs distribtution 

p(r) = Z- 1 e-^ = J2p S m\m){m\ (15) 

m 

after coupling to the heat bath for a time much longer 
than the relaxation time 7 -1 of the considered system. 
The occupation probability in the instantaneous eigen- 
state |m) is 

-p 3 E m (T 3 ) 

P S m(t) = Y' ' (16) 

which depends on the spectral structure and the tem- 
perature T 8 (s = h,l) of the relevant heat baths. Here 
(3 S = 1/(KT S ), and K is the Boltzmann constant. The 
partition function is defined by 

z* = Y, e ~ PsEm{Ts) ■ ( 17 ) 

m 

For the cyclic nature of QHE, the energy level spacings 
return at different instants. Then the net work done 
during a cycle can be calculated as 

AVF = Y, bm(n) - P l m (r S )] A£ roi (18) 

m 

where 

AE m = E m (T h )-E m {T l ). (19) 
Therefore, the PWC can be explicitly expressed as 

Aw = E [^( r i) - p L M] AE ™ > °- ( 2 °) 



V. CLASSIFICATION OF POSITIVE WORK 
CYCLES FOR 3-LEVEL QHE 

It is also known from Refs. 0, ^2 that the PWC for 
the 2-level QHE working between the source and sink at 
temperature and TJ can be written as 

A h 

^>T™ (21) 
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FIG. 3: (color on line) A sketch of an adiabatic evolution of 
the level spacings of a 3-level system during the step 2. A^ 
and A\(i — 1,2) are the two level spacings when coupling 
to the source and the sink. A h (A l ) is the maximum level 
spacing, i.e., A s = Af + A|, (s = h, I). 



Only when this PWC is satisfied can the positive work be 
extracted. This condition is counter-intuitively different 
from that T/j > 1} for a classical heat engine. Now, we 
naturally ask a question: Is this condition l(2*T|) universal 
for a multi- level QHE? In this section we will prove that, 
if the energy levels change properly, the PWC for a 3- 
level system can be looser than that for a 2-level system 

CD- " 

Our model considered here is a 3-level system with 
the adjustable level spacings Af and A 2 (s = h,l), as 
illustrated in Fig. 3. Here, we denote the ground state, 
the first excited state and the second excited state with 
subscripts 0, 1 and 2 respectively. We also introduce the 
dependent parameters A s = Af + A|, (s = h, I). For this 
kind of 3-level systems, if the PWC l(2Tj)l can be reduced 
to 



T h > Tt (A h /A l ) 
T h > 7] (Af/Ai) 0' 



(22) 



meanwhile both A 1 / A h < < 1 and A^/Aj- < 0' < 1 are 
satisfied, then we say the PWC is looser than that for a 2- 
level case since the PWCs for both the two substructures 
are looser than that for a 2-level system 1)21(1 . The two 
substructures are formed by combining level with level 
1 and level 2 respectively. We will prove in the following 
there indeed exist such cases for a 3-level system. This 
means that the PWC for a 3-level system can be improved 
in comparison with that for a 2-level case when the levels 
change properly. 

For the 3-level case we rewrite the expression l(18ll of 
AW in a more explicit form 



AW = 



G(A,A 1 ) + G(A 1 ,A) 



(23) 



Z h {T h )Z l {Ti) 

where Z S (T S ) is the partition function for the three level 



G(A, Ai) 
6(A,Ar) 



9(A,A 1 )(A /l -A'), 



(24) 

(25) 



Obviously AW is completely determined by T s , Af and 
A 2 (s = h,l), which is independent of the details of the 
changing of the level spacings. Thus, we only care about 
the initial level spacing A^ 1 and the final level spacing A' 
in considering the thermodynamic cycle. 

However, at arbitrary finite temperature, the PWC 
(|2*U|l is too complicated to be understood directly, so we 
switch to considering its high temperature limit. In this 
limit the PWC (|2*U1) can be simplified as 



L 



i-h A 1 



(A? - Ai) > 0, (26) 



where we have introduced 



F&ft) = [(2£-l) + (2 -£)»?]■ 
A' 



L = 



mZ h {T h )Z\Ti) 



in terms of an independent set of parameters 



^2 



A 1 

^2 . 



A^ 



(27) 



(28) 



A 



^4; A 



A\ 



According to the requirement of the quantum adiabatic 
evolution, the level spacings Af and A| are always 
kept positive during the two quantum adiabatic process. 
Therefore, the parameters r\ and A range in the interval 
(0, 1), i.e., < 77, A < 1. 

In principle, we can classify the multi-level QHE ac- 
cording to the changes of the level spacings in the ther- 
modynamic cycles. As to the 3-level QHE, there are al- 
together four sorts of operations corresponding to the 
following four cases: 



(I).A\ - A[ > 0, 
{II).A\ - A[ > 0, 



(III) .A'l 

(IV) .AI 



A! < 0, 
Ai < 0, 



A2 - A l 2 
A h 2 - A< 
A h 2 
A'i- 



> 
< 
- A!, > 
Ai < 



(29) 



as schematized in Fig. 4. 

Before analyzing the four cases in detail, we intuitively 
consider the physical mechanism that a 3-level QHE may 
better the work done. Namely, the PWC can be relaxed 
compared to that for a 2-level case. For instance, if the 
lower level spacing does not change during step 2 (see 
Fig. 5 a), A^ = A[ (dashed line), the 3-level QHE re- 
duced to a 2-level case, then the PWC cannot be bet- 
tered. However, when it comes to A\ > A[ (solid line), 
the lowing of energy level 1 produces extra work to better 
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4 



R(x,) R(t 2 ) R 

(III) 




R(x,) R(t 2 ) R 

(IV) 



Case (I), the above PWC (|30|l can be further simplified 
as 



T h >T l ^ r e. 

Here, the dimensionless parameter 9 is defined by 



(31) 



(32) 



Aiming to construct a 3- level QHE with a relaxed con- 
straints 11' I'll on temperatures, we need to analyze in what 
conditions 9 < 1 is satisfied. Obviously, as to 9 < 1, i.e., 
(2 — £) A < (2 — £) rj, there are two solutions. Solution I: 



£ > 2 and X > r/, 



(33) 



FIG. 4: (color online) Sketch map of four sorts of 3- level QHE. 
There are altogether four kinds of 3-level QHE, as illustrated 
above, according to the changes of the level spacings during 
step 2. Here (I), (II), (III) and (IV) are corresponding to case 
(I) to case (IV) in Eq. (29) respectively. 



and Solution II: 





R(x,) R(x 2 ) R 



FIG. 5: (color online) A 3-level QHE built up by adding a 
third level to the 2-level QHE. Intuitively speaking, it can 
make the PWC looser (better the work extraction) when the 
new level helps to produce positive work (a), but makes the 
PWC worse when it cost energy (6). 



the PWC. On the contrary, if the level spacings changes 
as that in Fig. 5 b, A\ < A[ ( solid line), the raising of 
the energy level 1 makes the PWC become worse. These 
two cases are corresponding to Case (I) and Case (III) 
in Fig. 4. We expect that the 3-level QHE of Case (I) 
can better the work extraction but Case (III) can not. 
We will try to prove this result and determine the phys- 
ical parameters that can better the work extraction in 
the following detailed analysis. Firstly we consider Case 

(I). 

Case (I): Physically, Af - A{ > and A£ 
mean both of the two level spacings decrease adiabati- 
cally during step 2. Since A^ - A^ > 0, the PWC 
can be simplified as 



A(, >0 



F^,r ] )T h >T l F{i,\) 



A 1 



(30) 



By noticing that £ > 1, F(£,rj) > 1 and F (£, A) > 1 in 



1 < £ < 2 and A < n. 



(34) 



For simplicity, we define a set of independent physical 
parameters in terms of the ratios of three level spacings 
to A'l, 



l _ Ai , 
Then it is easy to see that 



A', 



Ah 
^2 



(35) 



(36) 



are implied in our presupposition of Case (I). 

Based on the above analysis, Solution / can be explic- 
itly obtained to be 



to > r\ r. 



rf, > 1 



(37) 



We combine Eqs. I|36|) and (|37|) and plot them in a 
3-dimensional figure in coordinates of r[, r\ and r\. 
The solution is a 3-dimensional domain enveloped by a 
curved surface, r\ = r\r^ (denoted by A) and two planes 
r% = 1 — t[ + 7*2 (denoted by B) and r[ = 0. Any rep- 
resentative point in this domain can determine a 3-level 
QHE, which can extract positive work under the condi- 
tions T h > Ti (A h /A l ) 9 and A l /A h < 9 < 1. 

We further consider whether Solution I is consistent 
with the second inequality of Eq. I|22|) . Fortunately we 
can easily find such 9' smaller than unity for A\l A\ > 
A h /A l in Solution /. Thus, according to the the defini- 
tion (|22|l . it is proved that the PWC for 3-level system 
can be looser than that for a 2-level case if the level spac- 
ings change properly. 

Similarly we rewrite Solution 7T (|34|l in terms of r[ , r\ 
and rj, 



rU < r\ r. 



r 2 < 1 

2 „h 



(38) 
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FIG. 6: (color online) Diagram of solution I of case I. For a 
given Af , when A[, A2 and A§ are properly chosen such that 
the point (r[,r 2^2) is in the enveloped range, positive work 
can be extracted under a looser condition. 



FIG. 8: (color online) Sketch map of the net work AW as 
a function of Th for a given T. Part a illustrates a classical 
Carnot engine, while part b illustrates a 2-level QHE. 





r h 
'2 


1 1 








r\ 
2 \ 









(color online) Diag 



FIG. 7: 

representative points in the enclosed range are the results of 
inequalities l|38|l . and the physical meaning is similar to that 
of Fig. 6. 



Together with H36|) . we illustrate the result in Fig. 7. 
The physical meaning is the same as that of Solution I. 
But Solution 77 is not an ideal solution like Solution J, 
for we cannot find such solution that satisfy the second 
inequality. 

In order to make our result clear, we illustrate the 
above results in comparison with that for a classical heat 
engine and a 2-level QHE in Fig. 8 and Fig. 9. As to a 
classical heat engine, e.g., a Carnot engine, the PWC is 
Th > Ti, and the net work AW is a linear function of Th 
for a given T h But for a 2-level QHE, the PWC is Eq. 
(|21|l . and the net work is not a linear function of Th, as 
illustrated in Fig. 8. 

As to Solution I of case I of the 3-level QHE, the PWC 




a 





AW 












T 





T, T* / 
r 'A', / 


'A' 





b 



FIG. 9: (color on line) Similar to Fig. 8. Part a illustrates 
Solution I of Case / of the 3-level QHE, and Part b illustrates 
Solution II of Case I of the 3-level QHE. 



is looser than that for a 2-level case, because the PWC 
for both the two substructures are looser than that for 
the 2-level system, as illustrated in Fig. 9a. However, 
as to Solution II, the PWC is not looser than that for a 
2-level case, for the PWC for one substructure is looser 
than that for the 2-level system, but the other is not, as 
illustrated in Fig. 96. 

Case (II): In this case - A[ > and A' 2 l - A l 2 < 
0. Physically, this means that the upper level spacing 
increases while the lower level spacing decreases during 
step 2. We will prove in the following in this case there 
exist no such solution about Aj 1 and Af that can better 
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the work extraction compared to a 2-level case. 

To prove this conclusion we need to distinguish the 
following 4 situations: 



a 
b 
c 
d 



F(Z, V ) >0,F(£,A) >0, 

F(t, v )>o,F(s,\)<o, 
F(t, v )<o,F(s,\)<o, 
F(t, v )<o,F(t,\)>o. 



(39) 



a: The PWC can be reduced to Eq. fHJ|. Similar 
to above analysis, if we want to find a looser PWC, 9 
should be smaller than unity, i.e., 9 < 1. By noticing 
£ < 1 in Case (II), 9 < 1 can be further simplified as 
A < Tj, i.e., 



Ai _ A| 
A?T 



(40) 



However, this inequality is incompatible with the two 
constraints of Case (II): A\ - > and A£ - A^ < 0. 
Thus, in this situation, the PWC is no looser than that 
for a 2-level QHE for 9 > 1. 

b: The two inequalities F(£, 77) > and F(£, A) < 
can be reduced to 



I — 277 , 1 - 2A 



2-77 " 2 - A 
Thus we can obtain the inequality r\ > A, i.e., 



(41) 



A] 
A? 



A 



h > Alt' 



Similarly, this inequality is incompatible with Aj 1 — A^ > 
and A% - A l 2 < 0. Thus, in Case (II) F(£, rj) > 
and A) < can not be satisfied simultaneously. 
Therefore, there does not exist the PWC looser than that 
for a 2-level QHE in this situation. 

c: The PWC can now be reduced to be 



T h < T x 



A h 
AT' 



(42) 



In summary, we cannot find a PWC looser than that 
for a 2-level case in Case (II). After similar analysis 
about Case (III) and Case (IV), we find there is no 
desired solution in Case (III) and Case (IV) either. 
In conclusion, only in Case (I) can we find such desired 
solutions that the 3-level QHE can better the work ex- 
traction compared to the 2-level case. This conclusion 
agrees with the result obtained from our foregoing intu- 
itional consideration. 

We also would like to mention that under the criterion 
(^21 wc find the PWC for a 3-level QHE can be looser 
than that for a 2-level case. If we use other criterions, the 
result may be different. For example, if we use A2 instead 
of A in the first inequality of (|22|l . we can not find such 
3-level QHE that whose PWC is looser than that for a 
2-level case. Actually, for Case (I), the PWC (O can al- 
ways be simplified to inequality (J2TJ T h > T t (A h /A l ) 9. 
It can be proved that the coefficient (A h /A 1 ) 9 satisfies 



Af- Aff 
ApAl 



A h 
A*~ 



< max 



A| Aj 

A[ ' A l 2 



(44) 



This conclusion is also true for n- level QHE. Namely, if 
the PWC for a n-level QHE can be expressed as Th > Tk, 
then it can be proved that the coefficient k satisfies 



Ai 1 



A h 

71 

AI 



^ k ^ max 



A? 



A h 

n 

AI 



(45) 



Here A% and A' are the energy gaps between the ith 
and the (i — l)th energy level in Steps 1 and Step 3. 
Specifically, when all the level spacings change in the 
same ratio, i.e., A\/A\ = A%/A l 2 = ■•■ = A'^/A l n , no 
matter how the spectral structure is, the PWC for such a 
system has the same form as that for a 2-level case Th > 
Ti [A h / A 1 ) . The harmonic oscillator and the infinite well 
potential are two good examples [|| Q . 



VI. ILLUSTRATION: 3-LEVEL ATOM WITH 
DARK STATE 



By noting that we have assumed beforehand Th > T. 
Only when (A h /A l )9 > 1, the QHE is probable to ex- 
tract positive work. From Eq. (|39|l . F(£, rj) < is always 
satisfied in condition c. Then (A h /A l )9 > 1 can be re- 
duced to 



i = «-i>' 



!<»• 



(43) 



Obviously this inequality has no solution to the real pa- 
rameter £. Therefore the positive work can not be ex- 
tracted in this situation. 

d: The PWC can now be reduced to T h < Ti(A h /A l )9. 
The left hand side in the inequality is positive while the 
right hand side is negative (for 9 < 0). It is obvious that, 
there does not exist a PWC looser than that for a 2-level 
QHE in this situation, either. 



In the preceding section, we found in proper conditions 
the PWC for the 3-level QHE can be looser that that for 
a 2-levcl QHE (condition I). In this section, we will use 
a concrete example to demonstrate our results. We con- 
sider a toy model with 3-level system couples to a classical 
single-mode external field. The levels are adjusted by the 
external light field, which plays a similar role to an ideal 
piston in classical heat engine. In this ideal case we need 
not to consider the work for the controlling field to ad- 
just the level during the adiabatic process. Just through 
this external field the work (either positive or negative) 
done on the 3-level system means the change of energy 
of the 3-level system. Namely, in our theoretical studies, 
no extra work is required for the controlling field since 
we regarded the field as a part of the external entries . 
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The model Hamiltonian [16( reads 

H = HS \e) (e\ +n{\e) (1| + |e) (2| + h.c) , (46) 

where 6 is the common detuning (see Fig. 10). We have 
set the eigenenergy of the two degenerate ground states 
1) and |2) as zero, and that of the excited state |e) as 
huj e . £1 is the complex Rabi frequency associated with the 
coupling of the field mode of frequency lo to the atomic 
transition |e) — * |1) (|e) — * |2)). The detuning of lo is 6, 

LO = u> e - S. (47) 

We solve the eigen-equation, and obtain the eigenvalues 



E 
E+ 



(48) 



= o \ M + V ( hS > + 8 \ n 



- hS 



{hsy + 8 |n| 



Then this QHE has two level spacings: 

Ai = E - = i ^(M) 2 + 8|fi| 2 - m\ , (49) 

A2 = E + — En = 



The two level spacings are functions of two independent 
parameters S and |f2|. During Steps 2 the level spac- 
ings change adiabatically from A h (A J) to A' (A^) when 
the two parameters change slowly from (S h , \Q h \) to (S l , 
I ft 1 J ) . According to the systematical analysis in last sec- 
tion, this model can indeed serve as an improved 3-level 
QHE if the two parameters change properly. We plot the 
energy levels A (6, \Q\) and Ai (5, in Fig. 11. 

In a thermodynamic cycle, the level spacing changes 
adiabatically between Aj 1 and A' during the two adia- 
batic processes. For convenience, we define 



K s 



\9J 



(50) 




FIG. 10: (color online) The sketch of a 3-level QHE model. 
A A-type 3-level system interacting with a single classical ra- 
diation field. |e) is the exited state, while |1) and |2) are the 
two degenerate ground states. 



A. 2 




FIG. 11: (color online) The two level spacings as functions of 
hS and |fi|. The upper curved surface represents A s while the 
lower one represents Af. 



The two constraints of Case (I) in Eq. (|29H can also be 
expressed in terms of K s and 8 s : 



K l - M 



K h - HS 



k 1 + h8 l 
K h + hs h 



< 1. 



(51) 



Substituting Eq. (03 into Eq. (|57Jl . we get the Solution 
I in terms of K s and 5 s , 



8 l K h > S n K l , 
S h > 6 l . 



(52) 



Combine Eqs. I|51() and (|52|l . we find the inequalities hold 
only when <5 1 > 5 > 0. Thus the solution (|52|) can be 
reduced into a more compact form 



s h >s l > 0, 
n l /5 l < n h /s h 



(53) 



During step 2 of a thermodynamic cycle, we assume 
the initial level spacings are A h and Aj" corresponding 
to two positive parameters S h and \& h \- Then the level 
spacings change adiabatically to A' and A^ when the 
two parameters change slowly to (6 , |^'|)- The above 

result 1|53J) shows that, only when the final point (5 l , 
locates in the shaded triangle in Fig. 12 can we make the 
PWC looser compared to condition I for a 2-level case. 

Similarly, we substitute Eq. (03 into Eq. (|3^|) and 
then obtain Solution 7T in terms of K s and 8 s . 



S h <5 l < 0, 

n l /s l < n h /s h 



(54) 



The physical meaning is almost the same as that of Fig. 
12. For the initial parameters 6 and \Q h \, only when 
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FIG. 12: (color online) Illustration of Eq. J^HJ- In this so- 
lution the QHE can extract positive work while the PWC is 
looser comparing with that for 2-level QHE. 
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FIG. 13: (color online) Illustration of Eq. 1541 . The physical 
meaning is similar to that of Fig. 12. But this is not an ideal 
solution. 



the final point (S l , locates in the shaded triangle in 
Fig. 13 can the first inequality of (|22|l be satisfied, but 
the second is still not. Thus Solution II is not a desired 
solution. 

Before concluding this section, we would like to address 
a crucial problems in the physical implementation of the 
above "dark state model" . It is obvious that we do work 
on the system when varying the parameters S and in 
step 4. But how is the work extracted out of the system 



can not be imagined directly. Here, what we discussed 
in this section is only a toy mode to illustrate the basic 
spirit of physics in the quantum thermodynamic cycles. 
The work done by the system in step 2 can be indirectly 
understood as the consequence of the decrease of the sys- 
tem energy, which can be explicitly calculated from our 
foregoing analysis (see Eq. JJJJ) in section II. 

VII. CONCLUSION AND REMARKS 

We have quantum-mechanically formulated the work 
done and the heat transferred in the thermodynamic pro- 
cesses in association with the microscopic quantum tran- 
sitions. A class of QHE were universally proposed by 
using a multi-level quantum system as the working sub- 
stance and by deriving the thermodynamic adiabatic pro- 
cess from the quantum adiabatic process. We classified a 
3-level QHE based on the changes of the level spacings, 
and found when the parameters (and thus the level spac- 
ings) change properly, a 3-level QHE can better the work 
extraction compared to a 2-level case. 

Before concluding this paper we would like to point 
out that the two-parameter QHE proposed in the last 
section is only a toy model and one can not over esti- 
mate it. We have to say that the Hamiltonian for such 
an atomic system in the laboratory frame of reference 
is time-dependent and changes fast. Therefore one can 
roughly regarded it as something relevant to the sys- 
tem energy. It is still an open question to find a multi- 
level system with level spacings changing independently 
in practice. 

We also remark that the discussion about QHE is es- 
sentially semi-classical because we quantize neither the 
heat bath nor the controllable external fields. To build 
a totally-quantum theory for QHE, one need to use the 
generalized master equations with and without memo- 
ries. There are some interesting results in Refs. [l7l, fl8| . 
How to develop our present studies within this theoret- 
ical framework is to be considered in our forthcoming 
investigations. 
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